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1. INTRODUCTION AND SUMMARY 

Background 

Over the past several years , MIT ' s Laboratory for Information and 
Decision Systems (LIDS) has been conduotin*^ research for NASA on the 
properties of multivariable digital control systems [1] . These types 
of systems are becoming increasingly important as smalls powerful, flight- 
qualified digital computers take over the burden of control law implemen- 
tation in various NASA vehicles axid other control system applications. 
Examples include the shuttle orbiter, the HIMAT and P-8C DPBW aircraft, 
satellites and space probes such as Viking and Galileo, various proposed 
large space systems, and many more. 

The overall goal of the research program has been to evolve improved 
design methods for multivariable sample-data control laws. Research 
effort was concentrated initially on the primary available synthesis tool— 
namely the sample-data (discrete-time) Linear-Quadratic (LQ) regulator 
problem [Athans, 2]. Various properties of this problem formulation were 
studied, and key features of its solution were investigated. In the latter 
category, the basic "robustness properties" (e.g. multivariable gain 
margins and phase margins) of sample-data LQ solutions received particular 
emphasis. 

Two major conclusions emerged from these initial studies. Pirst, the 
sample-data LQ problem, like other sample-data synthesis formulations 
suffers from the over idealized time sequence representation of the control 
process. The process is daseribed only "at the sampling instant," with 
no information about intersample behavior or intersample control requirements 



other than what can bo predicted from values of states and controls at 
the sampling times [Levis, 31. The need for devices such as prefilters 
{which appear in virtually every practical digital control implementation) 
and the potential benefits of higher-order hold devices cannot be addressed 
within this pure sample-data framework. 

The second major finding concerns the robustness properties of 
sample-data LQ control systems. These turn out to be fundamentally inferior 
to their continuous- tijne counterparts. Recall that the latter enjoy impres- 
sive uncertainty tolerances including -6db to + «db gain margin or j^GOdeg 
phase margin in all combinations of control channels [Safonov and Athans, 4]. 
These margins are guaranteed for every continuous-time LQ design, inde- 
pendent of specific plant or cost function parameter values. In the sample- 
data case, the corresponding uncertainty margins are generally less than 
the ranges above and are functions of specific plant and cost function values. 
Hence, no a-priori robustness guarantees can be stated for sample-data LQ 
control laws [Safonov, 5] . 

Motivated by these two apparent limitations of the existing LQ 
synthesis methodology, the research effort was re-directed toward more 
fundamental issues of digitally implemented control systems. The first task 
of the redirected effort was to find a mathematical representation which 
properly captures both the continuous- t.lme (emalog) and the discrete-time 
(digital) processes which occur simultaneously in a digital control system. 

Such a representation was developed and is called the "hybrid operator model" 
of the control process. This model provides an analog input-routput yiew of the 
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control process which explicitly includes sampling operations, digital 
calculations, hold operations, and continuous plant evolutions. The 
structure of this operator is summarized briefly in Section 2 below and 
in more detail in a Master's thesis by A. Kostovetsky [6]. 

An immediate application of the hybrid operator is to explain the 
common use of prefilters in practical digital control systems. Simple norm 
calculations in [6] show that the hybrid operator has unbounded gain (in 
an appropriate function space sense) as the sampling process tends toward 
the ideal impulsive sampling normally assumed in sample-data theory. 
Physically, this means that it pro'vidss arbitrarily large amplification 
for certain inputs (e.g. noise). Non-impulsiva sampling, as obtained with 
pre-filters, bounds this amplification. 

The second task of the redirected research made use of the hybrid 
operator model to answer the following very basic approjcimation <juestioni 
How well can digitally-implemented control laws mimic analog ones? More 
specifically, if samplers, holds and digital algorithms are all selected 
to best approximate a given linear, time- invariant analog system, how good 
can the approximation be? The answer to this question is elegantly simple; 
The digitally- implemented system can exactly duplicate the impulse response 
matrix, G(t-Q), of the analog system at all points in the t, 6-plane 
except on a strip of width T (sample time) along the main diagonal (t«6) . 
Inside this strip, the hybrid system's impulse response must be zero on 
various triangular segments* Accordingly, this region of approximation 
has been named the "triangle strip." Details of the optimal sampling, hold, 
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and digital function for this approxiination result are again summarized 
in Section 2 and derived in detail in [6] . Its significance is that it 
provides a simple and clear picture of the basic limitations inherent in 
digitally- implemented controls* Such controls are fundamentally inferior 
to their analog counterparts because they cannot utilize all the input 
data in the triangle strip. As shovm below, this limits bandwidth, 
restricts performance, and precludes robustness guarantees such as those 
enjoyed by the continuous- time LQ regulator. 

New Research Results 

Motivated by the above characterization of hybrid system limitations, 
the research effort was continued to explore the qualitative and quantitatiye 
ways in which these limitations manifest themselves in hybrid system design. 
The results of this new research are reported here. 

The new research examines hybrid operator models from the viewpoint 
of "conic sectors theory" [Zames, 7, and Safonov, 5]. This viewpoint 
encompasses the approximation results above and provide an effective way to 
express their implied limitations in terms of conventional amalog control 
system concepts (e.g. bandwidths, loop gains, gain auid phase margins, etc.). 
The basic idea is to approximate hybrid operators (which may be viewed as 
periodically time-varying analog systems) by time-invariant linear analog 
models. These models form the "centers" of conic sectors, and their modelling 
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errors deteripine the sector '•radii''. The size of sector radii/ in turn, 
determine bandwidth limitations/ stability robustness reductions, and 
performance reductions, in accordance with existing sector-based control 
system design technic^^les iDoyle and Stein, 8, Lehtomaki, 13], 

The relevant conic sector conctipts and their use in analog system 
design are briefly developed in Section 3, They are then applied in 
Section 4 to hybrid operators, T'wo expressions are derived which define 
sectors radii as functions of the true hybrid operator and its time- 
invariant analog approximations. These expression show explicit depen- 
dences on the selected sampling rate, the selected sample- and hold filters, 
and the selected digital algorithm. These parameters may be chosen to 
minimize the radii (in an appropriate function space sense) and, thus, to 
minimize the bandwidth, robustness and performance limitations which the 
hybrid operator imposes. 

The conic sector expressions from Section 4 are illustrated with a 
small design example in Section 5. While this example are limited in 
scope it serve to illustrate the use of sector-based hybrid system ap- 
proximations and point the way to a general control system design philosophy 
which incorporates digital implementations under a more general common 
umbrella of feedback design for systems with approximation errors. 
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2. HYBRID OPEHATOR MODELS 


Ws will consider digitally-implemented control systems which can 
be represented by the block diagram of Figure 1. The three main func- 
tions associated with the controller block in this diagram are: 

1) The sampling operation which converts M-dimensional 
analog inputs u(t) on the interval (A-l)T<t££Tinto 
N-dimensional discrete samples Z"l,2, , , . , 

2) the digital algorithm which converts the N-dimensional 

sequences into L-dimensional sequences rij^/ 

k“0,l,. and 

3) the hold operation which converts the L-dimensional 

sequences into R-dimensional analog functions v(t) 

on the interval kT£t<_(k+l) T. 

The system’s sample time will be designated by the symbol t. These 
three functions will be assumed to have the forms 




f^(0)u(e)d9 » 


k 

I 

Z=1 




Zr 

IA-Dt 


f (e-^.T)u(0)d6 

o 


( 1 ) 

( 2 ) 


v(t) = = h^(t-kT)rij^ (3) 

The first of these equations is a simple analog convolution operation 
with weighting fxinction (impulse response) f^(X). This could be the 
weighting function of an analog prefilter, an approximate impulsive 
sample, or various other vector valued input averaging operations. 
Some examples are given in [6] . The second equation is a standard 







digital convolution with ooamoiant. Tha third i. a gan.ralitad 

output hold operation with weighting function This could he a 

simple constant to represent the common "aero-order hold»" hut in general 

it will he selected to achieve broader goals. Some examples are again 

given in [63 , Note that the controller is completely characterised by 

the two matrix-valued functions h (A) and by the coefficient 

0 0 

matrices 


Hybrid Operator Hepresentation 

Given the above description of a digitally-implemented controller, 
it is straight-forward ([6] , Section 2) to write its input-output 
operator representation, G, in the terms of an impulse response matrix, 
G(t,9), That is. 


V “ G u (4) 

*v 

where v and u denote functions on [0,«) related by the convolution 


with 


v(t) 



G{t,B)u(9)d0 


(S) 


k 

G(t,9) « h^(t-kT) I D. .f (0-31T) (6) 

Here k is understood to be the largest integer lesJH than or equal to 

t/T. We will refer to this input-output description of the controller as 

"the hybrid operator model" or simply as the "hybrid operator". Note that 

it is a time-varying linear dynamic system characterized by h , f ^ and 

o o 



I 


r-. 

»'i<i .Jt 
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Optimal Hybrid Approximation 

consider now the problem of finding a hybrid operator model 

6(h ,f »D, o) to approximate a continuous-time linear dynamic control 
0 o 

law with impulse response matrix 


G(t,e) 

where A,B, and C are given system matrices. 

Let the approximation criterion be to minimize 



(7) 


( 8 ) 


where v(t) and v(t) are the outputs of the hybrid and pure analog 
controllers, respectively, when excited by the same white noise input. 
Then it is shown in [6], Section 4, that the optimal approximating hybrid 
controller has the following sampling function; 


f (X) » 
o 

It's corresponding hold function is 
h^(X) * Ce^, 

and the digital algorithm is 



gAT (k-il) 


(9) 


( 10 ) 


( 11 ) 


Moreover, these parameters cause (6) to duplicate (7) exactly every- 
where except on the "triangle strip" of Figure 2. Note that the sampling 
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FIGURE 2 : Trianule Strip. 
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and hoXd function* (9) -CIO) of thi* optimaX hybrid approximation ara 
thamsaXvQs n-th ordar dynamic system*/ where n is the dimension o;' A« 

Henca/ the ovaraXX hybrid controXXar can be visuaXised as shown in 
Figure 3. 

As indicated earXisr/ the /iignificanca of the above rasuXt is not 
the OptimaX structure in Figrure 3 itseXf (after aXX/ the sampXing and hoXd 
functions are quite oompXax/ each XiteraXXy dupXicating the anaXog system ) , 
but rather the fact that the inherent hybrid system iimitations are so 
simpXy and cXearXy dispXayed by the triangXe strip in Figure 2. It foXXows 
from this figure that the minimium apjiroximation error is given by the 
error operator 

e (G-^)u (X2) 

vrhere has the impuXse response representation 

t 

e(t) ■ f 5(t-0)u(0)d0 (X3) 

QuaXitatively therefore / the hybrid system suffers an inherent time 
varying “data lapse*' with a maximum duration of T seconds (average x/2 ) , 
and with data weighting proportional to the desired impulse response » G. 
Hence, both the nominal function G and the sample time X contribute to 
the significance of the error. Small errors ara assured if G(l) is small 
over the whole interval 0£x<x and u(0) is relatively "smooth." These 
observations are given further interpretation below. 
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at t»T/ 

Sample Digital algorithm Hold 


FIGURE 3; The Optimal Hybrid Compensator. 




3. CONIC SECTOR CONCEPTS 


Tha kay faaturo of feha hybrid oporat:'?r modal 0 is that it rapraaonta 
an analog input-output viaw of digitnlly-iinplemantad aontrollarg. That is 
tha oparator maps continuous-time input fui^ctions, u {i,a, u(t), 0;<t<<») 

into continuous- time output functions, v (i.a. v(t), 0<t<«) , Ml intar- 

sample behavior such as limit cycling and aliasing is thus included in the 
raprosantation. Tha input-output viewpoint also it possible to 

utiliao certain function space notions of approximation — in particular, 
the concept of conic sectors. 


Conic Sectors t5,7] 

In order to discuss conic sectors, it is necessary to review a 

number of basic concepts from functional analysis. The first of these is 

the function space, i^. This is the collection of all n-dimansional 

vector- valued functions which are square integrable on 10,*”) . If the 

function x belongs to i”, than its norm is defined by 

2 



A 



1/2 

(t)x(t;dt < «> . 


Moreover, the inner product for any two functions x, , 
defined by ^ 


<x. ,x > 


/ =*1 


,^(t)X2(t)dt, 


(14) 

x„ in i” can be 


(IS) 


It follows that 
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The ohiaf limitation of the space for control system analysis 
is that it contains no unstable functions (i.e. fionctions with ijx||»^). 
This can be remedied by introducing the so-called extended space/ 

This space is the collection of all functions which are square integrable 
on all finite time intervals/ i.e. it contains all functions x which 
satisfy 


X , 

'T 


[/ 


x*(t)x(t)dt < « for all T 


(17) 


t 1 

Functions such as x(t)«e are included in , for example/ while func- 

2b 

tions such as x(t) « tan t are not. 

An operator such as G is a mapping which associates each functions in 
its domain {the set of input functions) with exactly one function in its 
range (the set of output functions). For our purposes/ the domain of G 
will be the space and its range will be some subset of We will 

also assume that G is causal and L^^-stable. Causality means that the 
output of G at time t»t^ does not depend on future inputs/ 

L -stability means that the norms of the input and output functions of G 
are related as follows j 


I |(5u| L < k I |u| I (18)' 

for soma k<<», all u in , and all T<» 

This implies that all bounded input functions produce bounded output func- 
tions. Note that the two norms on the right and left of this expression 

are different because they are defined on different 



OF POOR QUALITY 

lU J* 

function spaces (L^ and , Their ratio is often used to define yet 

2e 2e 

another norm for the operator itself, namely 




(19) 


This is the so-called induced operator norm, induced by norms on 
ana 1.^^. It ia a aorimon abuae of notation in the fnnotional analysis 
literature to use the same symbols, ] |*||, for function norms and for 
induced operator norms. The distinction is made clear by the arguments 
used with the symbol. 

Given these preliminaries, the operator G can be viewed as a subset 
of the crossproduct space . This is shown schematically in 

Figure 4, where input functions are ’ represented by points along the x-axis, 
and the operator is represented by the graph itself. Within this pictorial 
framework, a conic sector is a cone-shaped subset of x which 

contains many potential operator graphs. In mathematical terms., it is the 
collection of all operators, H, which satisfy 

<Hu-{C+R)u, Hu-(C-R)u> < 0 , for all T<« and all u 6 

*^26 

( 20 ) 

Here <•,•> is the inner product on L and C and R are causal, 

X 26 

L_ -stable operators. C and R define the boundaries of the conic sector 
26 • 

and are called the "cone center" and "cone radius", respectively, for obvious 


geometric reasons. 



17 - 


original PASS IS 
OF POOR QUALITY 



FIGURE ii Conla Sector (C,R) 
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Using elementary meuiipulations of (20) , it is possible to show that 
conic sectors provide a way to approximate G by C. This is done as 
follows: Suppose (20) is satisfied by H"G [i.e. G is a member of Sector 

(G,R)]. Then 

<Gu “(C+R)u, Gu -(C-R)u> 

M Arf M jjl 

» <(G-C)u - Ru, (G-C)u + Ru> 

« <(G-C)u, (G-E)u> - <Ru,Ru> <0 

for all T<« and all u e (21) 

2 © 

Hence, 

I I (G-C)ul I < I iRul L (22) 

•I ~ ' T ■" - T 

for all T«» and all u 6 L”!" 

r, 2e 

This shows that the output of C approximates the output of G within 
approximation error bounded by the norm of the radius operator R. 

Note from (22) that sector membership also implies that 

I I (G-C)V\ 111 |RP| I (23) 

for any operator V, 

Sector-Based Feedback Design 

In the last several years, the conic sector concept has been 
recognized as an important tool in feedback analysis and design 
[Safonov, 5, Zames,7, Doyle and Stein, 8, Lehtomaki,13] . The basic idea 
is that very complicated plant operators G (perhaps nonlinear, infinite 



-X9“ 


dimensional, time-varying, etc.) can be reliably approximated by simple 
cone centers C (usually linear time- invariant finite dimensional systems) , 
provided that the approximation error is properly accounted for in the 
design process. This "proper accounting" usually weans that a design 
based on C must be restricted to maintain stability robustness, and its 
performance goals must be relaxed to account for inherent degradations. 

These restrictions generally increase as the magnitudes of approximation 
errors grosy. The basic relationships between design restrictions and conic 
sector approximation errors are summarized below. 

stability Robustness Restrictions — We will treat the generic feedback 

problem where G is the loop siTransfer operators of a unity feedback loop, 

as shown in Figure 5. G is assumed to be any causal, L 2 ^~stable member of 

Conic Sector (C,R) for which the feedback loop is "well-posed" [Willems, 9). 

-1 

This means that (I+G) exists and is causal. We will also assume that 
the nominal feedback operator (I+C) ^ is -stable. Then the loop's errors 
in response to disturbances, d, and reference commands, c, are given by 
A 

e = G - y 

fV Arf 

« (I+G)"^(c-d) 

« (l+C+{G-C)]“^(c-d) 

•s* 

=» (I+C)“^[I+(G-C) (I+C)”^]-^(c-d) 


(24) 


20 
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*. Y 

^ € Sector (Ct/^) 


FIGURE 5; Generic Feedback Loop. 
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All inverses in this expression exist and are causal by assumption. 

-1 

Moreover, (I+C) is stable by assumption. It then follows that e will 
be bounded and the feedback loop will be stable in the sense of (18) if 
and only if the operator 

[I+(6-C) (I+C)"^]"^ 

is stable. As shown in [Sandell, 10], this stability requixwent is 
assured whenever 

I I (G-C) (I+C)“^| I < 1 (25) 

or, using (23) , whenever 

I lR(I+C)“^| I < 1 (26) 

or equivalently, whenever 

I 1rG"^[C(I-^C)“^3 1 1 < 1 (27) 

and C ^ exists. 

Equation (27) shows that stability can be achieved for all plants 
in the Conic Sector (C,R) if the feedback system's nominal closed loop 
responses, C(I+C) are restricted to be small for all inputs which have 
large normalized conic sector approximation errors, RC A frequency 
domain interpretation of this restriction is given shortly. 

Performance Degradation -- Looking beyond mere stability, equation (24) 
can also be used to show that sector approximation errors impose inherent 
degradations in the performance of feedback systems. In accordance with 



our input-output viewpoint, we wiXl evaluate performance in terms of 
the error magnitudes generated for some specified subset of disturbance 

y» 

and/or command f motions. That is, let S C ^ specified set of 

functions (say, all sinewaves with frequency less than w ) and let q>0 

o 

be the smallest scalar such that 


I Ifl It " II I I^ 1 'll Is-jl It 

for all (c-d)e S and all T<<» 

A# 


(28) 


Then q will be taken as a quality measure for the feedback system' a 
performance. 

We will assiiine that the nominal system has quality measure 
with q^ designed to be sufficiently small. The question then is to 
determine the actual value of q which applies in the presence of sector 
approximation errors. While it is generally difficult to compute this 
value exactly, equation (24) can be readily manipulated to obtain a 
useful upper bound. In particular, from (24) 


e ^ (I+G) (c-d) + (I+C) Ac 


(29) 


with 

Ac = {[l-XG-0 (I+C)“^3"^-I>(c-d) 

« -{[I+(G-C)(I+C)"^3"^(G-C)(I+C)“^}(c-cl) (30) 

Af fV 

Since we have already established conditions for the inverse in the last 
expression to be stable, it follows that under those same conditions 



1 

1-1 |R(ItC)”^' 


||R(I+C)"Nc-d)||^ 


(31) 



and using (31) in (29) m hava 


ll;llTi, ||(i+C)'^(o-a) 


Ilin-cr’-|| 

1-||R(I+C)"*I1 


iiR(i+cr^{o-d)iij 


(32) 


This equation requires careful interpretation. Note that the first term 
on the right hand side is simply the performance of the nominal system. 
Hence, for inputs restricted to the specified set, (c-d)€ S, we have 

A# Ak* 


||(H-C)"^(c-d)||^ < qjlc-all^ • <33) 

The second term on the right hand side has three elements! 

A 

Element (1)! | j (I+C)“ 1 1 

This is the (global) norm of the nominal 

sensitivity operator and will be denoted by s* 

Note that s* can never be less than unity with 
strictly proper plants since for these there 
exist functions d such that Cd is arbitrarily 
small. In well designed systems, however, s* 
approaches unity quite closely. 

Element (2): 1-| |R(I+C)’”^| | 

This is a (global) stability margin of the feedback 
system and will be denoted by m*. From (26) , its value 
is guaranteed to be positive if the system is stable. 

Typical values for well designed systems hover around 

0.5. 

— 1 

Element (3): j jR(I+C) I ^T 

This is the function norm of the sector radius operator, 
operating on the nominal feedback errors. For inputs 



r. «i 
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f > « 




. \ !:' I 


reatrictad to lot the 3:adlus operator aatisfy 

I|R;II i'o ll^llT »“) 

for all Q e(I+C) ^ S and all T<«. Then 

A* 

I |R{I+C)'^=-a) 1 1 j I |c-d| (35) 

with r^^Q necessarily less than unity from (26) 

Putting these elements together allows us to reduce (32) to the following 
simple forms 

l|on^<a+ -^VqJ|o-d||j (36) 

for all (c-d)e S and all T<«> 

M A# 

'This show that sector approximation errors degrade feedback performance 
in a direct and simple way. The nominal quality measure is simply scaled 
upward by the factor ^^‘r^/m*. T)ia latter increases with increasing 
sector approximation errors on the function set S, and also with deter- 
iorating global sensitivity and stability margin properties of the nominal 
design. 

Frequency Domain Interpretations — If the cone parameters C and R are 
stadjle linear, time-invariant operators, then the norm inequalities in 
(26), (27) and (36) can be readily evaluated with classical frequency-domain 
methods. This is done with the following norm identity for linear time 
invariant operators [Desoer and Vidyasagar, 11] 

I|C|| » max <7[C(s) ] 
u 


(37) 
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C(s) is ths (matvi:<-vaiued} Laplaca transform of tha impulsa 
rasponsa, C(t), of tha oparator G, and a['J# £l'J danota its maximum 
and minimum singular values . Using (37 ) , it follows that (26) and 
(27) are satisfied if 


o(R(jw)3< £ti+C(jw)3 for all 0<o)«» 


o(Rc"^(jw)3< U(c(jw) (l+c(jw))”^] 
■« a[i+c“^(jw)3 


for all o<o)«» 


Note that (27)' imposes explicit magnitude constraints on the nominal 

frequency response. At all frequencies where the normalized error, 

-1 

RC (jw) / is large (compared with unity) , the inverse loop transfer matrix 
-1 

C (j«) must also be large, and hence, the loop transfer matrix itself 
must be small, since RC"‘^(jw) typically grows at higher frequencies, this 
constraint imposes explicit limitations on achievable crossover frequencies 
and on closed loop bandwidths. 

Note that conic sector approximation errors play identical roles in 

\ 

the above sense to any other modelling errors already recognized in the 
design process. 


The performance degradation equation (36) can be similarly interpreted. 


We get 


'We will use the same symbol (C) to denote both the impulse response and 
its transfer functions. The differentiation will be made clear by arguments. 



t 
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When substituted into (36), these equations show that performance 
degradation due to sector approximation errors ean be minimised by 
designing good nominal sensitivity and stability margin properties over 
the entire frequency range (small s*, large m*) and by assuming small 
sector approximation errors over w<a)Q (small r^) . 
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4. CONIC SECTORS FOR THE HYBRID OPERATOR 

Motivated by the above discussion of conic sector concepts in feedback 
design r this section proceeds to derive two expressions for radii of conic 
sector approximations of the hybrid operator developed in Section 2, The 
first expression is based on time-domain analysis and yields a conservative 
non-dynamic radius. The second expression is based on frequency-domain 
analysis and yields a sharper dynamic radius operator. Some examples using 
the second radius are given in Section 5. 


A Non-Dynamic Radius 

We found in Section 2 that the hybrid operator is a linear, time- 
varuing convolution operator with impulse response matrix G(t,0) defined 
by equation (6) . Our objective now is to approximate this operator by a 
linear time- invariant convolution with impulse response C(t-6). The ap- 
proximation error should be such that it can be bounded by a sector radius 
operator as in equation (22) . 

To achieve the latter objective, consider the convolution operator 
for the eri^or, i.e. 


with 


e = (G-C)u 

t 

e(t) = J* [G(t, 0)-C(t-0) ]'u(0)d6 
0 


= / AG(t,0)u(0)d0 , 

Q 


t>0 


(42) 
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* t 


This expression can be bounded as follows? 

T 


I G-C « su] 


sup 

-ii4i, 


T, l|u||j»l 


L*'o 


(t) e(t)dt 


T 

t 

t 

y* dt 

Jae 

JdX u'^(e)AG^(t,8)AG(t,X)u(X) 


0 0 


sup J* 


L T t t 

fit fie f 

0 0 0 


dX |u^(6)|ai;i\G(t,e)]a[AG(t,X;3 |u(X) 


(43) 


Now let m(t) be a scalar-valued function such that 


a[AQ(t,e)] < m(t-6) 


vt,e 


The right hand side of (43) is then bounded by 


t t t 


" T L 0 0 0 


dX w(0)m(t-Q)m(t-X)o)(X) 


Here M is a scalar convolution operator defined by 
V - M(ji) 

with 


/ 


v(t) =■ / m(t-0)w(6)d9 , t>0 

'0 


(44) 


(45) 


(46) 
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Using (37) / it now follows that a bound on the error operator (42) can 
be expressed in terms of the Laplace transform, m(s), of (46). That is, 


1 i^“^l 1 £ 

Sr 

U7) 

w 

s=j(ja 


Moreover, since m(t)^ 0 for 

all t, it can be 

shown that the maximum 

in (47) occurs at w=«0 and hence 

00 


r = m(s) = I 

f m(t)d0 

(48) 

s=0 

r 

) 



In light of (47) , the operator 
V = Ru 

defined by 

v(t) = ru(t), t>0 (49) 

forms a valid radius for the sector approximation. 

Note however, that the calculation of this radius if formidable. 
It requires the following steps; 

1) Impulse response evaluation 

AG(t,6) = G(t,6) - C(t-0) 

2) Singular values calculations 


a[AG(t,0)] 


Vt,0 
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3) Function maximization 


0[AG(t,6) ]< m(t-6) 


Vt,0 


4) Integration 


/” 


r « j m(t)dt 
0 


The most difficult of these steps is 3) which requires that we find an 
integrable bounding function for <7 [AG] over the entire t, 0-plane. This 
step is made substantially easier by noting that the error operator is 
periodic along the (t+0)-axis whenever the digital convolution coefficients 
.0^^ in equation (6) satisfy \ We then need to maximize 

only over a single period for each value of the time index t-0, as 
illustrated in Figure 6. 


A Dynamic Radius 

In addition to its computational difficulties, the radius operator 
(49) suffers from conservatism. Note that its impulse response is a delta 
function, and hence, its Laplace transform is constant. It therefore 
assigns the same conic sectors error levels to low frequency inputs as it 
does to high frequency ones. This is. inconsistent with experience. We 

. 7f 

know that a digital controller can approximate signals at w« “ with 

TT 

much greater fidelity than it can approximate signals at ^ • This 

emperical fact is not captured by (49) . 




f 




m(X) 


max AG(t,0) 
0_<vi<X 


repeats on each 
similarly shaded 
square 

t=(y+X)/2 

e=(y-X )/2 


FIGURE 6; Calculation of m(t) . 



An alternate radius which captures this frequency dependence has 
been developed by means of the following classical arguments. We start 
with the fact that the hybrid operator can be represented in the classical 
sample data, block diagram form shown in Figure 7. [Franklin and Powell, 12 J . 




D(z) 

A A A JkA A 

H(s) 

^ 

■ F(s) 

T 



Prefilter 

Sampler 

Digital 


Hold 



Computer 


- analog signals 

/vvs/w digital sequences 

FIGURE 7; The Hybrid Compensator • s Sample Data 
Block Diagram. 


While we cannot find a Laplace or Fourier transform for this operator, 
it is possible to find the transforms of the output fvmctions produced by 

particular inputs. Namely 


v(jw) = Y 


Lk““<» 


P(joj-jo) k)u(jtO“jaj k) 
s s 


(50) 


ICOT 

where D(e ) is the z-transfoinn of the digital convolution sequence 

evaluated at z=e^^"*', and where the infinite sum with to = ~ 

s T 

represents the effect of sampling. The analog signal before the sampler has the 
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Fourier transform F(jw)u(ja)) . After the sampler^ the signal is 
mathematically represented by a sequence of impulses whose transform is 
the infinite sum. For notational convenience, we will represent this 
sum simply by ( ^ F u ), where F. and u are F(jo)-jw k) and u(jta-jw k) , 

K K K K S S 

respectively. 

Using (50) in Parseval's theorem and assvuning that all functions are 

in L^, we can now write expressions for the hybrid operator's errors; 

2 


with 


(G-C)u 


i2 1 


2tT 


2^ 


00 


L 

CO 


v(j( 0 ) - c(jw)u(jw) 


£ I? “ 


dw 


dlt) 




“ hdf, 
T k 


r HOP -c 
T 0 


kf^o 


k=0 


Using Schwartz’s inequality, (51) can be bounded from above. 
2 


(G-C)u 


2TT 


/(iIgJ'Vi I-vIV- 

-CO yk 


(51) 


(52) 


(53) 
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rt-n 

. 1 , 




Now adding and subtracting HOF^l in the first sum of (53) and 
interchanging the order of integration and summation converts the bound 
to 


i- 

2tt 



\1 

't 


HDPq-cI^-17 HDFo 



(S4) 


In each of these integrals we now perform a change of variables, 

03«W' + 0 ) Jl. Note that this is a frequency shift by integral numbers 
s 

of sample frequencies, q . Hence, the periodic functions 

S 

0 0 

1d| and ( I |f. I ) remain unchanged. We also re-interchange the order 
k ^ 

of summation and integration. This gives 



27r 


2TT 


l( 


t Vh-h 


u 


do)' 


00 

filH 


loo 


D 


2 

.U dW 


(55) 


Note that the right hand side of (55) has the form 


2tt 


J|R|^|u|^dOl, 


which according to Parseval' s theorem, is equal to the norm Ru 
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Comparing this to our sector definition in equation (22) (restricted 
to functions) then shows that (55) represents a valid sector radius 
for our hybrid operator. The center of this sectors can be arbitrary, 
However, it is evident from (55) that the choice C* with Fourier 
transform 

C*(jw) « ~ H(jW)D(jw)F(jW) (56) 

generates the smallest radius, R*, with Fourier transform 



"I I I J |H(joJ+j03 )l)|^|D(e^‘^'^)|^lF(ja)+ja3 k)|^ (57) 


The sector defined by these expressions is explored by a simple example 


in the next section. 
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5. AN EXAMPLE 

Vhe second sector developed In Section 4 is illustrated here with 
a classical lead- lag compensator example, While this is a very simple 
example/ it captures the essential issues of sector-based approximations 
and serves to illustrate their potential utility. 


The Analog System 

We start with the analog compensator 


G^(S) 


lOCa+.ll 

s+1 


(58) 


We show how to implement this compensator using a digital computer, 
and then how to compute the center and radius of its Sector (C,R), 

The digital computer is embedded in the hybrid compensator as in 
Figure 7, We arbitrarily choose a sample interval of t » .6283 seconds, 
which places the sampling frequency at 2'iT/T « 10 radians per second, 
which is ten times the frequency of the pole location. The z-transform » 
D(z) (used for the digital computer) is computed by the classical pole- 

zero mapping technique tl2] . The pole at s=-l maps to a pole at 

ST ST 

z = 6 = .5335, and the zero at S!=»-.l maps to a zero and z?=e » .9391. 

The constant is chosen so that D(z) at z=l is equal to G (s) at s«0. The 

a 

result is 


D(z) = 7.660 


Z-.9391 


Z-. 5335 


(59) 
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The prefJLlter F(s) is chosen to he a low pass filter with a single pole 
at S"-tr/T» The result is 


P(s) 


S 

s+5 


(60) 


The hold device is chosen to be a zero-order-hold^ which has the 
Fourier transform 

i nm - . 


(61) 


Two Cone Centers 

According to equation (56), the cone center C*(jw) is the given by 


C* (jo)) 



e^‘^'^-,9391 5 

e^^''-.5335j 


(62) 


Note that this choice of cone center has the special property of being 
the "straight through path," formed by replacing the sampler in Figure 7 
with the gain term ~ as shown in Figure 8, 



FIGURE 8; The Straight Through Path. 
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Unfortonately/ this choics also has the undesirable property of being 
infinite dimensional and difficult to work with in state space realizations. 
For this reason, we will also examine an alternate center which approximates 
(62) by a low order state space model. This alternate choice is 



and is obtained by approximating the zero order hold with a first order 
filter and replacing D(e^^”^) with 


Cone Radii 

Following our previous derivation, the cone radius corresponding to 
center C*(s) is given, by equation (57). The radius for C(e), on the 
other hand, requires a slight modification of this equation to account 
for difference between c and c*. The appropriate modification is obtained 
from (55) and has the form 

Both equations (57) and (64) can be evaluated numerically by computing 
the infinite sums over k and H. It can be shown that these sums converge 
for all functions H, F, and C-C* which roll off as least as fast as 1/oj 
at high frequencies. In fact, if the roll-off is 1/w exactly, then the 
convergence properties of each sum are proportional to 


I ~ - 


Asi r 


6 


1.644934 


(65) 
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Wq hava emoiirlcally detorminsd that thase infinite sums can ba 
tarminated at plui and ndnuB 20 terms. Their general convergence pro- 
perties, however, still need to be better understood. 

Graphical Results 

The net results of all these calculations are best viewed as 

conventional Bode plots of the various functions involved. We start 

with Figures 9a and 9b, which show gain and phase characteristics of 

X 

the zero /^rder hold, ~ H, and its first order approximation. Note 

that the approximation is good up to about 2.4 rps, or roughly 

Figures 10a and b show Bode plots of the analog lead compensator, 

G (s) and the two centers, C(s) and C*(s). Note that these centers have 
a 

more phase lag than the analog compensator, due to the extra lag added by 
the prefiltor and holdi- However, C and C* again track closely up to 
about 2.5 r/s. 

The corresponding radii were computed according to equations (57) 
and (64) . One of the required infinite sums is shown in Figure 11. We 
see that above 2'ir/T * 10 radians the infinite sum is approximately unity. 
The radii themselves are plotted in Figure 12. We see that they agree 
closely for high frequencies (OJ^l rad/sec) but differ substantially at 
lower frequencies. This difference is caused by the terms in 

(64) . Note that the sum of these terms alias high frequency differences 
between C* and C down to lower frequencies. In particular, the differences 
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at aJ**w ■ 10 r/s appears as a DC level in the radius corresponding to 
C(s). This illustrates that any design advantages offered by a simple 
cone centers must be paid for with large radii and with the associated 
more severe st^lbility robustness restrictions and/or performance degra- 
dations discussed in Section 3. 

Some Basic Tradeoffs 

Whether or not a given sector approximation is useful for feedback 
design depends on the relative magnitude of its radius when compared 
with the corresponding center. We have already observed this relation- 
ship in Section 3. In particular, if a loop transfer operator G is any 
member of Sector (C/R), then our stability robustness condition (26) 
required that 

||R(I+C)‘^c-d)||^ <|l£-d||^ (66) 

for all inputs (c-d)e l!? and all T<<» 

«s# ^ 20 

Moreover, to get good performance we required that (I+G) ^ be small 
for some specified subset of inputs, i.e. 

l|(I+C)-^c-d)|l^ < ||c-d|4 (67) 

for (c-d)e S and all T<«® 

If q «!/ it follows that 

(I+C)"^(c-d) « C"^(c-d) (68) 

<s/ <v* 
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and substituting this into (66) we get 


llRC-N=-d)|4<||=-d||^ (69) 

for all c-d 6 S, all T<«®. 

*y »v 

This last expression shows that the sector approximation is useful 
for feedback design only if its normalized error magnitudes are less 
than unity for the signal set S. This turns out to be specific 
statement of a very basic truth in feedback control. We cannot close 
feedback loops with substantial gain unless the (normalized) loop 
uncertainties are less than one. 

In the example above, of course, we have computed radii and centers 
not for a complete loop transfer operator but for the hybrid lead com- 
pensator alone, h complete loop transfer operator using this compensator 
will be given by 


G 


loop 


- G P 


= [C+(G-C)]P 


(70) 


for some plant P and some compensator center C. The normalized loop 
transfer radius condition (69) then becomes 


|R, C'^ (c-d) 

I loop loop ~ ~ 


'T 


|RP(CP)-^(c-d)||^ 
|RC-^c-d)||^ <||=-d 


(71) 


for all (c-d)e S, T<« 
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This shows that a sector radius for the compensator alone is also 
useful for feedback design only if its normalized magnitude is lass 
than unity for the input set S, 

As in Section 3, Condition (71) can be interpreted in the frequency 
domain via 

atR(jw)C(jw)"^]<l (71)' 

for soma frequency rauage w£w . 


The Sample Rate Tradeoff — Condition (71) ' was evaluated for our 

second sector approximation of the lead compensator, as defined by 

(63) and (64) . The evaluations were made with several sampling times 

ranging from t=0.628 sec (w =10 r/s) to t» 0.0628 sec (100 r/s) . Results 

s 

are plotted in Figure 13. They show that Condition (71)' is not 
satisfied at any frequency for our baseline value T»0.628, However, 
it is satisfied over increasingly larger frequency ranges as t decreases; 


T Sampling Rate Frequency Range 


.628 sec 
.314 
.125 
.0628 


10 r/sec 
20 
50 
100 


none 

0.03 - 5. r/sec 
0 - 13. 

0 - 30. 


This result illustrates the dramatic effect that sample rate selection 
has on digital compensator design. Oiir hybrid implementation of G (s) , 
as modelled by Sector (C,R), does not yield a useful compensator unless 
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T is less than approximately 0.2 seconds. We note that this conclusion 
is quite consistent with classical "rule of thumb" for selecting sample 
rates. An old rule due to Stein, for example, is 

"... 5 to 10 samples per radian of tne highest 
significant frequency in the control loop,.," 

In our example, the highest significant frequency is perhaps 1.0 r/s 
(where the centers CSC*' stop developing load) . The resulting rule- 
of-thiimb sample time would therefore be 0.1-0. 2 sec. 

The Discretization Tradeoff — In addition to sample time, a second major 

design parameter in digital implementation is the discretization technique 

used to convert G (s) to D(z). In our results so far, we have used an 
a 

established pole-zero matching technique as our baseline. Figure 14 
compares normalized radii obtained with this technique against the following 
alternates, all at t«0.628(12]; 

1) the "forward rectangular rule" (replace s in 

G (s) by (z-l)/T) 
a 

2) the "backward rectangular rule" (replace s in 
G (s) by (z-1)/zt) 

A 

3) the "Tustin rule" with pre-warping about 
[replace s by (to^/tan w^T/2) (z-1)/ (z+1) 3 , 

The differences between these techniques are seen to be small. None produce 

a useful sector approximation for the compensator at this sampling time. 
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We coroplote our diacusaion of the lead compenaator example at thia 
point. It ia evident that the example and its two aimpl® tradeoff! 
above only aoratoh the surface of the potential utility and range of 
application which conic sector concepts offers for digitally implemented 
control system design. It is hoped that these results serve to 
motivate the needed additional research and design studies. 


-SS"* 

6, COMCLUSION 

This report has prsssntscl research resuXts obtained under NAS;v 
Grant No, NAGl-2j "Hybrid Operator Models £or Digitally-Implemented 
Control Systems." The results establish a new method of analysis for 
digitally-implemented (hybrid) control systems based on conic sector 
concepts from functional analysis. Conic sectors are used to ap- 
proximate the true complex time-varying nature of a hybrid compensators 
with a simple time-invariant analog system. The errors of this ap- 
proximation are rigorously accounted for in the analysis, both in terms 
of their effect on stability robustness and on performance. 

Two specific conic sector approximation were developed for general 
hybrid systems. One sector produces a conservative non-dynamic radius, 
while the other is less conservative and frequency dependent. The latter 
sector has so far been derived only for stable scalar systems. It was 
illustrated with a simple lead compensator example. This example confirms 
the ccmputational feasibility of conic sector analysis and also serves 
to illustrate its potential utility and range of application. 

The conic sector concept promises to unify both analog and digital 
control design techniques under the more general common umbrella of 
feedback design for systems with approximation errors. Much work remains 
however to achieve this end. Tighter, less conservative conic sector 
approximations for general multi-input multi-output hybrid systems remain 
to be discovered, multi-rate systems remain to be analyzed, and the entire 
area of efficient numerical algorithms for conic sector analysis remain to 
be explored. 
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